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We study the consequences of an electron-phonon (e-ph) interaction that is strongly peaked in the 
forward scattering (q = 0) direction in a two-dimensional superconductor using Migdal-Eliashberg 
theory. We find that strong forward scattering results in an enhanced Tc that is linearly proportional 
to the strength of the dimensionless e-ph coupling constant Am in the weak coupling limit. This 
interaction also produces distinct replica bands in the single-particle spectral function, similar to 
those observed in recent angle-resolved photoemission experiments on FeSe monolayers on SrTiOa 
and BaTiOa substrates. By comparing our model to photoemission experiments, we infer an e- 
ph coupling strength that can provide a significant portion of the observed high Tc in these systems. 
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Introduction — A flurry of scientific activities has been 
generated by the discovery of an enhanced superconduc¬ 
tivity in FeSe monolayers grown on SrTiOs (STO) sub¬ 
strates [1-20]. On its own, bulk FeSe has a modest super¬ 
conducting transition temperature Tc ^ 9 K [21]; how¬ 
ever, when a monolayer is grown on an STO substrate, 
Tc is increased dramatically [1]. Most reported Tc values 
cluster within 55 - 75 K, close to the boiling point of 
liquid nitrogen (77 K). (A surprisingly high Tc ^ 107 K 
has also been reported in in situ transport measurements 
[9].) This discovery has opened a pathway to high-Tc su¬ 
perconductivity through interface engineering, which has 
already produced high-Tc’s in systems such as FeSe on 
BaTiOs (BTO) [8] and FeTei_^Se^ on STO [22]. 

Determining the origin of the Tc enhancement in these 
interface systems is critical. At the moment, proposals 
include charge transfer between the substrate and FeSe 
[2-4, 20], electric field [6] and strain effects due to the 
substrate [5, 8], and phononic related effects such as en¬ 
hanced e-ph coupling in the FeSe layer [1, 13, 16] or across 
the interface [7, 19]. Strong evidence for the latter has 
been provided by a recent angle-resolved photoemission 
spectroscopy (ARPES) study [7], which observed replica 
bands in the single-particle spectral function of the FeSe 
monolayer. These replicas are interpreted as being pro¬ 
duced by coupling between the FeSe 3d electrons and 
an optical oxygen phonon branch in the STO substrate. 
Moreover, the replica bands are complete copies of the 
corresponding main bands, which implies that the re¬ 
sponsible e-ph interaction is strongly peaked in the for¬ 
ward scattering direction (small momentum transfers). 
Such momentum dependence is notable because it can 
enhance superconductivity in most pairing channels [23- 
31]. As such, this cross-interface coupling provides at the 
same time a suitable mechanism for the Tc enhancement 
in the FeSe/STO and FeSe/BTO systems [7, 8]. 


In this Letter we explore this possibility and exam¬ 
ine the consequences of strong forward scattering in the 
e-ph interaction for superconductivity and the spectral 
properties of a two-dimensional system. By solving the 
momentum dependent Eliashberg equations, we show 
that a pronounced forward scattering results in a Tc that 
scales linearly with the dimensionless e-ph coupling con¬ 
stant (see below) in the weak coupling limit. This 
is in stark contrast to the usual exponential dependence 
predicted by BCS theory. Furthermore, this coupling 
produces distinct replica structures in the spectral func¬ 
tion similar to those observed experimentally. By com¬ 
paring our model to experiments [7], we infer a signif¬ 
icant e-ph contribution to the total Tc observed in the 
FeSe/STO system with a modest value of A^. 

Formalism — To model the FeSe monolayer we con¬ 
sider a single-band model for the FeSe electron pockets, 
which includes coupling to an oxygen phonon branch in 
the STO substrate. The Hamiltonian is given by 

k,cr q 

+ Y 5(k,q)cbq_,Ck,,(6Lq + 6q), 

^ k,q,cr 

where cj^ ^ (<^kcr) ^q (^q) electron and phonon 

creation (annihilation) operators, respectively, is the 
band dispersion, (2^ is the phonon dispersion, and ^(k, q) 
is the momentum-dependent e-ph coupling constant. 

We calculate the single-particle self-energy due to the 
e-ph interaction using Migdal-Eliashberg theory. Us¬ 
ing the Nambu notation with fermionic Matsubara fre¬ 
quencies Un = (2n + 1)7t/(3, where P = 1/T is the in¬ 
verse temperature, the self-energy is Y(k, icj^) = ico’^[l — 
Z(k, icc;n)]To + x(k, ico’n)T 3 + 0(k, ico’n)Ti, where fi are the 
Pauli matrices, Z(k, icj^) and x(k, icj^) renormalize the 
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single-particle mass and band dispersion, respectively, 
and 0 (k, is the anomalous self-energy, which is zero 
in the normal state. In Migdal-Eliashberg theory, the 
self-energy is computed by self-consistently evaluating 
the one-loop diagram and is given by 

S(k,zw„) = X 

^ q,m 

f 3 G(k q, iuj 

m )^3 

where Z)(‘^)(q, iw^) = — is the bare phonon propa- 

gator, and G“^(k, iujn) = — E(k, iujn) is the 

dressed electron propagator. 

In what follows we parameterize the electronic dis¬ 
persion as = — 2 t[cos(/Ca^a) + cos{kya)] — /i with t = 
75 meV and /i = —235 meV. This choice in parame¬ 
ters produces at T an electron-like Fermi pocket with 
kp = 0.97/a and a Fermi velocity vp = 0.12 eV-a/h 
along the ky = 0 line, where a is the in-plane lattice 
constant. This closely resembles the electron pocket at 
M point measured by ARPES experiment. Since first 
principles calculations indicate that the relevant oxygen 
phonon branch in STO is relatively dispersionless near 
the T-point [14, 32, 33], we approximate the phonon with 
a flat Einstein mode = Q = 100 meV {h = 1), which 
is consistent with the observed energy separation of the 
replica bands [7]. Furthermore, as we are interested in 
the case of forward scattering, we neglect any potential 
fermion momentum dependence in the e-ph interaction 
and set ^(q) = exp(—|q|/go), as microscopically de¬ 
rived before [7, 19]. Here, qo sets the range of the cou¬ 
pling in momentum space. For different values of qq we 
adjust go to obtain the desired value of the dimension¬ 
less e-ph coupling constant A^, which is computed from 
the Fermi surface averaged mass enhancement in the nor¬ 
mal state Xm = (— L-q) [^^]' Throughout we 

assume an s-wave symmetry for the gap function, con¬ 
sistent with the observations of a fully gapped state on 
the Fermi level [1, 7, 10, 18]. Finally, we neglect the 
Coulomb pseudopotential /r*. One can therefore regard 
our Tc values as upper bounds for the e-ph contribution 
to the FeSe/STO system. 

Analytical Results — Before proceeding to full numer¬ 
ical solutions, we can gain some insight by first consid¬ 
ering the case of perfect forward scattering, where the 
e-ph matrix element is a delta function \g{q)\‘^ = g^S^N 
with = Xm^‘^ [34]. In the weak coupling limit, we 
further set Z(k, = 1 , x(k, = 0 , and therefore 
(/)(k,= A(k,With these approximations, the 
gap function on the Fermi surface is given by 



FIG. 1: (color online) The superconducting gap at the small¬ 
est Matsubara frequency A(z7r//3) as a function of tempera¬ 
ture for various values of the e-ph coupling strength Am, as 
indicated. The e-ph coupling constant ^(q) is strongly peaked 
in the forward scattering direction with = 0.1/a. The inset 
shows Tc as a function of Am, which is extracted from the data 
in the main panel. The thin dashed line is the result in the 
limit of perfect forward scattering (see text). The shaded area 
represents the values of Am that are relevant for FeSe/STO 
[34]. 


tion is linearlized by setting Aq = 0 for T ^ Tc and we 
set Un=i/^ = 0. This results in the condition for Tc 

XrfjiTl^ \ 2^4 

Pc ^ + cj^) ■ 

The Matsubara sum can be performed exactly, yielding 
our final expression 

^ ^mPc 217 + ^4 cosh(^4/dc) — {^1 Pc) sinh(^4/3c) 

^ ^ 1 + cosh(f^/3c) ' 

For the case of FeSe, Tc and the hyperbolic func¬ 
tions dominate. To the leading order, the critical tem¬ 
perature is quasi-linear in the coupling strength in the 
weak coupling limit, Tc = 2 + 3 ^ similar result was 

obtained in Ref. 23 in the context of the cuprates us¬ 
ing square-well models.) For Xm = 0.16 and ^4 = 100 
meV one obtains Tc = 75 K, which is a remarkably high 
temperature for such a modest value of A^,- 

The increased Tc should be compared to the standard 
BCS value obtained for a momentum-independent cou¬ 
pling. In this case, the linearized gap equation simplifies 
to [34-36] 


1 = irTcXm T, 

\uJm I <^^ 1 ) 



— Xn 


In 


A 

27rTcJ 



A / . \ _ Xm^ ^{i^m) 2^4 

To determine Tc we take the ansatz /X{iujn) = Ao/[l + 
(cjn/f^)^] and follow the usual steps [36]: the gap equa¬ 


where we have expanded at large Qp/Tc and 'ip{z) is the 
digamma function [34]. This form produces the usual 
exponential behavior for the critical temperature, Tc = 
l.lSftp exp{—l/Xm), which predicts d, Tc = 2.5 K for 
A^ = 0.16 and ^4^ = 100 meV. 
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FIG. 2: (color online) The temperature dependence of the spectral function for several values of the e-ph coupling Am- 


Comparing these two results, one sees that the origin 
of the enhanced Tc lies in the momentum decoupling [24] 
that occurs in the Eliashberg equations when the inter¬ 
action is strongly peaked at q = 0. In the BCS case, the 
integration over the Fermi surface is equally weighted at 
all momenta, leading to a term in the BCS gap 

equation and subsequently a leading logarithmic behav¬ 
ior. In the forward scattering case, there is no integra¬ 
tion over momentum so the term remains, resulting 
in a leading behavior that scales like 1/Tc [34]. Thus, 
strong forward scattering serves as an ideal mechanism 
for producing high-Tc superconductivity [30]. Further¬ 
more, a strong forward scattering peak in the coupling 
constant means that this interaction will contribute in 
most pairing channels [7, 23-31]. It can therefore act in 
conjunction with other active unconventional channels, 
providing another means to increase Tc further. 

Numerical Results for Tc and the superconducting gap 
— In real materials the e-ph interaction is expected to 
have a finite range go in momentum space [7, 33]. There¬ 
fore we now consider an interaction with a finite width 
by numerically solving the full Eliashberg equations for 
an e-ph coupling constant g(q) = go exp(—|q|/go). Fig. 1 
shows the superconducting gap at the lowest Matsubara 
frequency A(ki?, Att//^) as a function of temperature for 
several values of and go = 0.1/a. We find that the 
superconducting Tc is already large for a modest value 
of \m and increases approximately linearly with A^, in 
the weak coupling limit; however, the finite range of the 
coupling in momentum space reduces the total Tc slightly 
with respect to the perfect forward scattering limit (see 
the inset of Fig. 1). The linear dependence of Tc with 
respect to A^ may account for the wide variation of re¬ 
ported Tc values in the literature, as differences in sample 


preparation or doping are likely to result in differences in 
the screening of the e-ph coupling and subsequently Tc. 

Replica Bands — The above results show that, in prin¬ 
ciple, a modest coupling to a phonon with a forward 
scattering peak is capable of accounting for the large Tc 
enhancement observed in FeSe on STO and BTO. The 
natural question is then how much of the experimental 
Tc is accounted for by this coupling? The observed shape 
and intensity of the replica bands [7, 8] provide us with 
a direct means to estimate this by comparing our model 
to experiment. To do so, we calculate the single parti¬ 
cle spectral function M(k, ca) = —ImGii(k, Cc;)/7r, which 
requires the analytic continuation of the self-energy to 
the real frequency axis using the method of Ref. 39 (see 
also [34]). Fig. 2 plots the temperature evolution of the 
spectral function obtained from a full numerical solution 
to our model for several values of A^, as indicated on the 
left, and go = 0.1/a. In all cases clear replica bands are 
produced by the coupling, offset in energy from the main 
band by a fixed energy, which is Q for small values of A^. 
The separation, however, grows for increasing A^. This is 
due to x(k, cj), which shifts the main band upward in en¬ 
ergy. This is most clearly seen in the A^ = 0.33 results, 
where the value of kp has visibly shrunk in the main 
band. In addition, for stronger values of A^ we begin 
to see the formation of a second replica band located at 
^ 2Q below the main band. Thus the observation of only 
a single replica band in the bandstructure of FeSe/STO 
is consistent with a small A^. 

An intuitive picture for the intensity and energy split¬ 
ting of the replica band can again be obtained in the limit 
of perfect forward scattering. On the real axis, the zero- 
temperature self-consistent equation for the self-energy in 
the normal state can be written as = goG(co’ + ^4). 
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FIG. 3: (color online) The spectral function for a momentum 
at the band bottom (k = 0 in our model, the M point in the 
experiment) for T = 30 K, = 0.1/a and Am = 0.02 — 0.22. 
The key feature of the forward scattering mechanism is the 
appearance of the mirror band {Z-) next to the main band 
(Z+). The relative separation Au and intensity = 

A(0, Ct;_)/A(0, Ct;+) of these two features is shown in the inset, 
and increases approximately linearly with Am- The dashed 
lines show the corresponding result in the perfect forward 
scattering limit and the shaded area represents the values of 
Am that are relevant for FeSe/STO [34]. 


For ^ 0“, the lowest-order solution is 
[note that the Ck 7^ 0 solution can be obtained by shifting 
the self-energy S(k,u;) = S(u; — ^k)]- The poles of the 
Green’s function are at cj = S(u;), which has the solution 
uj± = — § ± + 4^q. The spectral weight of each 

pole is given by Z± = [1-|^U=u;±]“^ = [1+ 

For small A^ = average energy 

separation between the poles is Auj = f2[l+2Am + 0(A^)] 
and the ratio of the spectral weight is ;^ = Am + 0(A^), 
thus providing a direct measure of A^- 

The spectral weight ratio and energy splitting between 
the main and replica bands can be extracted from our nu¬ 
merical simulations for finite values of go- Fig. 3 shows 
A(k, uj) for k = (0, 0) as a function of A^ with go = 0.1/a. 
The behavior matches our expectations gained from the 
perfect forward scattering limit: both the distance be¬ 
tween the bands and the relative spectral weight grow 
with increasing A^,, though the rate of increase is slower 
than for the case of perfect forward scattering. ARPES 
experiments on the FeSe/STO system [7] observe a spec¬ 
tral weight ratio of ^ 0.15 — 0.2 [34]. Comparing to our 
model calculations, we extract a value of \m ^ 0.15 — 0.2. 
This corresponds to a Tc ^ 60 — 70 K and a gap mag¬ 
nitude of A ^10 — 15 meV, which are consistent with 
measurements [1, 7, 10, 18]. 

In Fig. 4 we present the evolution of the spectral func¬ 
tion for increasing values of go where \m is fixed to give 
the same value of Z_/Z+. As expected, the replica bands 
are observed to smear both in energy and momentum as 
the value of go is increased. This shows that a weakly 


go = 0.1/a go = 0.3/a go = 2/a 
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FIG. 4: (color online) The spectral density A(k,ct;) along 
the k = (^/a, 0) cut for go = 0.1 (left), 0.3 (middle), and 
2 (right). In all three panels the temperature is T = 30 K 
and Am = 0.14, 0.125, and 0.25 in the left, middle, and right 
panels, respectively. 


momentum-dependent coupling (large go) to an optical 
mode does not reproduce the observation of a perfect 
replica band, with the same effective mass and termina¬ 
tion points in the Brillouin zone. Consequently, strong 
forward scattering is a necessary ingredient to under¬ 
stand the experimental observations [7]. 

Summary and Conclusions — We have examined the 
consequences of e-ph coupling that is strongly peaked in 
the forward scattering direction on the spectral proper¬ 
ties and superconducting transition of a two-dimensional 
electronic system. We demonstrated that such a cou¬ 
pling produces distinct replica bands in the electronic 
bandstructure consistent with recent ARPES measure¬ 
ments on EeSe/STO and EeSe/BTO interface systems. 
In order to reproduce the experimentally observed spec¬ 
tral function, we find that relatively modest values of the 
e-ph coupling are needed with A^ ^ 0.15 — 0.2. Strong 
forward scattering results in a momentum decoupling of 
the Eliashberg equations, which subsequently produces 
a larger superconducting Tc in comparison to the pre¬ 
dictions of conventional BCS theory. As a result, the 
inferred values of A^, predict Tc values on the order of 
60 — 70 K due to e-ph coupling alone. 

We stress that our results do not exclude the presence 
of another unconventional pairing channel such as spin 
fluctuations. The predicted values of Tc and A will be 
reduced somewhat by the inclusion of the Coulomb pseu¬ 
dopotential /i*. This reduction, however, can be over¬ 
come by the combination of the e-ph and unconventional 
interactions, since forward scattering will contribute to 
Cooper pairing in most channels [7]. An obvious way to 
distinguish between these possible scenarios is to measure 
the oxygen isotope effect. If a purely phononic mecha¬ 
nism is present then Tc should have an isotope coeffi¬ 
cient a = —91og(Tc)/91og(M) = 1/2, while the energy 
separation between the replica bands should decrease by 
^ 0.5(18 — 16)/16 ^ 6% for rich substrates. Alter¬ 
natively, in a multi-channel scenario, the isotope coeffi¬ 
cient a will be reduced from 1/2 when the unconventional 
channel is significant in comparison to the e-ph interac- 
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tion [37]. This provides a clear means to distinguish be¬ 
tween these scenarios. 

Finally, we note that e-ph coupling with a pronounced 
forward scattering peak has been studied in several con¬ 
texts related to of unconventional superconductivity in 
the cuprates [24-30] and pnictides [31]. Moreover, it is 
also now being addressed in the context of nematic fluc¬ 
tuations [40, 41]. This suggests forward scattering has a 
broader applicability in enhancing superconducting be¬ 
yond interface systems. 
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ELIASHBERG EQUATIONS 

In this work we have solved the momentum-dependent Eliashberg equations and obtained the electron Green’s 
function on the real axis, following the method in Ref. 1. Here, we briefly outline the main aspects of this approach. 

In the main text we introduced the electron self-energy in the Eliashberg equation on the imaginary axis with the 
Nambu notation 


t{k,iuJn) = iuJn[l - Z{k,iuJn)]ro + x(k,icJn)T3 + (t){k,iuJn)ri, (1) 

where fi are the usual Pauli matrices, Z{k, iuJn) and x(k, iuj^) renormalize the single-particle mass and band dispersion, 
respectively, and 0(k, iuJn) is the anomalous self-energy, which is zero in the normal state. In Migdal-Eliashberg theory, 
the self-energy is computed by self-consistently evaluating the one-loop diagram as shown in Eig. 1 and is given by 

S(k,iw„) = -ia;„)f3G(k + q,za;„)f3, (2) 

^ q,m 

where is the bare phonon propagator, G~^(k,iuJn) = — E(k, is the dressed 

electron propagator, N is number of momentum grid points, and [3 = l/{k^T) is the inverse temperature. Eor 
the dispersionless optical phonon mode we study here, Qq = Q. Eor forward scattering the e-ph matrix element is 
independent of the fermion momentum k and is given by 

|^(k,q)p = |^(|q|)P=^oV(|q|), (3) 

where /(|q|) is a function peaked at |q| = 0. Eor the perfect forward scattering, we take /(|q|) = YJq where Jq is the 
Kronecker delta function. Eor a general forward scattering we take /(|q|) oc exp(— 2|q|/go) with a proper normalization 
factor so that the result in the go ^ 0 limit can be compared with the result obtained using the Kronecker delta 
function. 

To obtain the electron self-energy for real frequencies, one needs either to perform an analytical continuation by 
Fade approximants, or to solve the Eliashberg equations directly on the real axis. The latter method is cumbersome 
while the former is unreliable; however. Ref. I introduced a method whereby the real-frequency self-energy can be 


S(k, iUn) 



FIG. I: The electron self-energy within Midgal-Eliashberg theory is given by the one-loop diagram above, where the Green’s 
function in the diagram is not the usual bare one but the fully dressed Green’s function. By solving the self-energy self- 
consistently, we implicitly include higher order terms in the one-loop diagram, but the vertex corrections are neglected according 
to Migdal’s theorem. 
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found iteratively by using the results on the imaginary axis. Given the imaginary axis Green’s function G(k, 
the real frequency self-energy is found by solving the iterative equation 

A/i X 1 fi I /I \ 1 ^3G(k + q, 

+ V ^ l5(k, q) |^B(z)f 3 G(k + q,uj-z + ir])f 3 ^ ^tanh + coth , (4) 

where B{z) = 6{z — ft) — 6{z ^ ft) is the spectral function of the phonon, and > 0 is a small number. 

For a specific momentum dependence /(|q|) in the coupling function |^(k, q)p = ^o/d^D? coefficient go 

to fix the total dimensionless coupling constant A^, formally defined as (see, for example. Ref. 2), 


Ak = —Re 


aE(k,a;) 


— (-^k) — 


9^ 1^=0 

I]k<^(^k)Ak 

Sk'^(^k) 


= Re[2’(k,w) 


(5) 

(6) 


For the relevant low temperature range Tc ^ T ^ ^2, the temperature dependencies of the self-energy in the normal 
state and A^ are negligible, so we calculate A^ at T = 100 K, if not otherwise specified. Here, (...) denotes the 
average over the Fermi surface, and ^(^q) is approximated by a Gaussian 


<5(3:) = 21" . (7) 

Xm is a measure of the mass enhancement due to electron-phonon coupling, m*/m = 1 -h A^, where m* is the effective 
electron mass. 

For a normal metal (that is, with a smooth coupling function ^(q) and ft Eq- (4)-(6) gives the textbook 

result for the total electron-phonon coupling constant A 

Am « A ^ E l5(k - k0|25(ek)(5(a'), (8) 


where N-p = is the density of states per spin. For the forward scattering case, however, the mass 

enhancement Xm following Eq. (6) differs significantly from the A defined in Eq. (8) for a normal metal. Eor the 
forward scattering with a width go in momentum space, Xm ^ QoX as go ^ For one set of parameters used in our 
numerical simulation, g = lb meV, N = 64x64 discrete momentum points and go = 0.1/a, we obtain A^ ^ 0.24A from 
numerical calculation. The Eq. (8) is therefore unsuitable to characterize the strength of electron-phonon coupling 
for forward scattering and to predict the corresponding Tc due to electron-phonon interaction. Consequentially, we 
use Xm as the strength of electron-phonon coupling in the main text. 

Numerically, we solve the Eliashberg equations with a momentum grid of = 64 x 64 grid points, a cut-off energy 
of 600 meV for the Matsubara and real frequencies with the real frequency discretized into 1 meV intervals, and a 
small-scale parameter = 15 meV. Eor a given set of parameters, we start the computation with initial self-energies 
Z = 1, X = 0 and a small anomalous self-energy <p = 0.7 meV. The Eliashberg equations are then iterated, where 
upon each iteration the chemical potential is adjusted as to fix the total number of electrons. When the difference 
between two consecutive self-energies is less than 10“^ meV, we have reached the convergence. 


BCS GAP EQUATION 


Eor comparison, we derive the standard BCS gap equation from the Eliashberg theory. An extended review of 
Eliashberg theory and the BCS limit can be found, for example, in Ref. 3. 

In the weak coupling limit, we can neglect the normal electron self-energy by setting Z(k, iuj^) = 1 and x(k, iuj^i) = 0. 
The anomalous self-energy is now equal to the gap function, 0(k, icjn) = A(k, For an isotropic 5-wave gap, the 
gap function is momentum-independent. Under these assumptions the Eliashberg equation reduces to 


q,m 


■^2 

?k+q 




(9) 
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Unlike for the forward scattering, we now assume that |^(k, — iuOrn) is momentum independent. Now 

the only momentum-dependent term in the Eliashberg equation is the electron dispersion ^k+q- In this case, the sum 
over momenta can be replaced by an integral over the free electron density of states A/'(e), 


-Y _ I _ 



1 


( 10 ) 


Since most of the relevant physics happens close to the Fermi level, we replace the N{e) by its approximate value 
A^(0). The resulting integration over electron energies changes the l/cj^ behavior into a l/lcj^l. 



1 


7rA^(0) 


1 

+A2(iw„) 


(11) 


The coupling strength as a function of energy can be expressed as = —N{fi)\g \^so that we recover 

the BCS gap equation. 


A{iL0n) = 'kT'Y “ i^^m) 


A{iWm) 


+A2(iw„) 


( 12 ) 


The critical temperature can be found by taking the limit = 0 and assuming that \{iuJn — and A(iuj^) are 
constants A and Aq with a hard cut-off at the Debye frequency The linearized gap equation with Aq = 0 now 
takes the form 


1 = ttTA Y 

l^m I 



A 




(13) 


where we have expanded at large and 'ip{z) is the digamma function. The critical temperature is given by the 

exponential function of the total coupling constant A 

Tc = (14) 


This exponential suppression can be traced back to the integration over the electron energies. This, in turn, is possible 
because the vertex and phonon part of the gap equation are momentum independent. When the vertex has strong 
momentum-dependence, the integration over electron energies cannot be performed and the leading behavior is still 
instead of l/lcj^l, as will be shown when we consider perfect forward scattering. 


PERFECT FORWARD SCATTERING 


For the forward scattering with ^(q) oc exp(— |q|/go), the Eliashberg equations can be solved analytically in the 
limit of perfect forward scattering, go ^ 0. The results of this limit are discussed in the main text; here we provide a 
more detailed derivation of the proper coupling strength A^ and the corresponding critical temperature Tc. 

The coupling for the perfect forward scattering is 

lfi'(q)P =fi'0^<^q ( 15 ) 

on a discrete lattice with total N grid points, where is the Kronecker delta function normalized to unity in the 
sum A = 1. (In the large N and continuum limit, A^q ^ (27r)^^(q), and the Dirac delta function ^(q) can 

be simulated with a peak function with finite width go? as for the general forward scattering case.) The one-loop 
electron self-energy now reads 

= — [ dz -- Ttanh-% + coth , (16) 

where B{z) = ^( 2 ; —U)—^(z + U). The subscript “(1)” of the self-energy means we have performed a single iteration of 
the solving the Eliashberg equation, where we substitute in the bare electron Green’s function instead of the dressed 
Green’s function. At the Fermi surface, the electron energy vanishes (^k = 0) so the self-energy simplifies to 


UJ 


ReE^^^(kF,co’) = gl coth 


n 


’ 


(17) 
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which is momentum independent. In the low temperature limit, coth(^4/2T) ^ 1, and by taking the derivative of the 
self-energy at cj = 0 we obtain the total coupling constant as defined in Eq. (6) 

The numerical value of Xm from the self-consistently calculated self-energy is very close to the value given by the 
above equation. For a nearly perfect forward scattering with a peak width go = 0.1/a, even after considering the 
small temperature dependence, the numerical value differs from the above result by a few percent. 

We now continue to derive the linear dependence of transition temperature Tc in terms of the coupling strength 
in the weak coupling limit for the perfect forward scattering. In this limit, we can neglect the normal electron 
self-energies, that is Z(k, = 1 and x(k, = 0. We assume an isotropic 5-wave gap. After summing over the 
momentum q, we obtain the self-consistent equation for the gap function 0(k, = A(k, 



A{tUJn) = 2^ - 


A(iuJm) 


2Q 


P ^ + /S?{iuJrn) + (^n “ ' 

LOm 

To find the critical temperature, we take an ansatz for the gap function with the form 

^0 




(19) 


( 20 ) 


Our numerical results verify that this shape properly represents the Matsubara frequency dependence of the gap. 
Setting uJn=i/Cl = 0 and linearizing the gap equation with Aq = 0, we find the equation for the critical temperature. 


Pc ^ w^[l + (w„,/f2)2](112 

UJm 


where the sum over Matsubara frequencies can be performed exactly by the method of contour integration. We finally 
obtain 


X^2Q -\-Q cosh ^ — 3Tc sinh ^ 
2Tc 1 + cosh 7p- 

^ c 


( 22 ) 


It is reasonable to assume that the critical temperature Tc is much less than the optical phonon frequency, as Q is 
typically of the order of 100 meV ~ 1000 K. In this limit, the hyperbolic functions dominate and cancel each other in 
the numerator and denominator of the above equation. We find a critical temperature of 


Tr = 


An 


2 T 3An 


-n. 


(23) 


This gives a remarkably high critical temperature, even for a modest coupling constant A^ because of the quasi-linear 
dependence. For example, Xm = 0.16 with Q = 100 meV yields a critical temperature of Tc = 75 K. 


SPECTRAL WEIGHT RATIO AND A^ FROM ARPES DATA 

By fitting ARPES data from Lee et al [4], we found that the ratio of spectral weights at the peak of the main band 
and the mirror band is approximately ^ 0.15-0.2, as shown in Fig. 2. This implies a value of the coupling 

constant Xm r\j 0.15-0.2. 
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FIG. 2: ARPES data (integrated and background subtracted) from Lee et al. [4] is fitted using the sum of two Lorentzians 
with height ratio 0.25, 0.2 and 0.15. The ratios in the range ~ 0.15-0.2 ht the experimental data best. The peak at 

— 100 meV energy in ARPES data is not included in our single band model calculations. 







